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THEORY  OF  THE  CURRENT-DRIVEN  ION  CYCLOTRON 
INSTABILITY  IN  THE  BOTTOMSIDE  IONOSPHERE 


1 .  Introduction 

Electrostatic  ion  cyclotron  (EIC)  waves  are  one  of  the  more  important 
plasma  wave  modes  in  the  near  earth  space  plasma  environment.  Recently, 
EIC  waves  have  been  invoked  to  explain,  among  other  effects,  strong 
perpendicular  ion  heating  leading  to  the  formation  of  ion  conics  [Okuda  and 
Ashour-Abdalla,  1981].  Ion  cyclotron  wave  phenomena  have  been  observed 
both  at  high  (several  Rg)  and  low  (<  1000  km)  altitudes  along  auroral  field 
lines.  Kintner  et  al .  [1973,  1979]  have  reported  H+  EIC  wave  observations 
using  S3-3  satellite  data  at  high  altitudes  in  the  auroral  zone.  At  low 
altitudes,  0+  EIC  waves  have  been  observed  [Ogawa  et  al.,  1981]  on  a 
sounding  rocket  equatorward  of  a  discrete  auroral  arc.  Yau  et  al.  [1983], 
also  using  sounding  rocket  data,  have  reported  particle  and  wave 
observations  of  low  altitude  (400-600  km)  ionospheric  ion  acceleration 
events  consistent  with  ion  cyclotron  wave  phenomena.  Bering  [1984]  has 
given  indirect,  but  compelling,  evidence  of  short  wavelength  (k  pi  >  1 
where  2ir/k^  is  the  density  fluctuation  scale  size  perpendicular  to  the 
magnetic  field  and  p  is  the  ion  gyroradius)  ion  cyclotron  wave  activity 
at  low  altitudes  (<  350  km)  in  the  diffuse  aurora  using  sounding  rocket 
data.  Evidence  for  short  wavelength  (k  pi  >  1)  ion  cyclotron  wave 
emissions  in  the  auroral  E-region  has  been  presented  by  Fejer  et  al,  [1984] 
using  backscatter  radar  observations.  Photometric  observations  in  the 
bright  auroras  at  -  130  km  altitude  by  Martelli  et  al.  (1971)  indicated 
excitation  of  EIC  waves  at  frequencies  lying  between  25  and  32  Hz, 
correspoding  to  gyrofrequencies  of  N0+  and  C^"*".  Bythrow  et  al.  (1984), 
using  Hilat  satellite  data,  have  concluded  that  the  very  large  earthward 
directed  currents  seen  near  the  equatorward  edge  of  the  diffuse  aurora 
exceed  the  threshold  for  the  excitation  of  ion  cyclotron  waves. 

Kindel  and  Kennel  [1971]  first  discussed  the  ion  cyclotron  instability 
in  a  space  plasma  context  by  showing  that  the  electrostatic  ion  cyclotron 
instability  has  the  lowest  threshold  among  various  current-driven 
instabilities  in  the  high  latitude  auroral  space  plasma.  They  only  briefly 
discussed  the  role  of  weak  collisicnal  effects  on  the  evolution  of  the  EIC 
wave  instability.  They  did  note,  however,  in  the  context  of  the  EIC 
Manuscript  approved  July  9,  1985. 
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instability  that  electron  collisions  were  destabilizing  and  ion  collisions 

were  stabilizing.  Physically,  the  ion  cyclotron  instability  is  related  to 

the  dissipative  effect  on  electrons  in  their  parallel  motion  which  impedes 

the  "instantaneous"  redistribution  of  the  electron  fluid  to  a  Boltzman-like 

distribution  in  the  wave  potential.  In  the  collisionless  limit,  the 

dissipation  is  due  to  the  wave-particle  (Landau)  resonance,  while  in  the 

collisional  case,  it  is  caused  by  the  electron  collisions  with  ions  (or 

neutrals).  The  results  of  Klndel  and  Kennel  [1971]  are  only  applicable  to 

the  collisionless  and  weakly-collisional  topside  ionosphere  [Chaturvedi  et 

al . ,  19841  .  Recently  Ganguli  et  al.  (1985)  have  proposed  a  mechanism  that 

depends  on  the  inhomogeneity  of  the  transverse  electric  fields  to  generate 

ion  cyclotron  waves  in  a  collisionless  plasma  in  the  absence  of  field- 

aligned  currents.  However,  several  of  the  aforementioned  experimental 

studies  have  indicated  ion  cyclotron  wave  emissions  can  occur  at  low 

altitudes  in  the  collisional  bottomside  ionosphere.  D'Angelo  [1973],  using 

the  work  of  Varma  and  Bhadra  [1964]  and  Levine  and  Kuckes  [1966],  has  shown 

that  the  current-driven  EIC  instability  might  be  excited  at  low  altitudes 

(E-region)  where  the  ion-neutral  collision  frequency  is  comparable  to  the 

ion  cyclotron  frequency.  This  work  has  been  extended  to  the  collisional 

auroral  F-region  altitudes  by  Chaturvedi  [1976]  based  on  the  analyses  of 

/ 

Chaturvedi  and  Kaw  [1975]  and  Mllic  [1972].  However,  Chaturvedi  [1976] 
considers  collisional  effects  on  the  ion  cyclotron  instability  using  a 
fluid  approximation  ( k^p ^  <  l).  This  analysis  is  not  applicable  to  the 
recent  experimental  observations  which  indicate  k  >  1 .  For  example, 
Beri_ng(  1984)  and  Fejer  et  al.  (1984)  have  reported  observations  of  meter 
size  irregular ities  corresponding  to  k^p^  *  1 •  Indeed,  the  transition 
between  collisionless  and  collisional  excitation  of  the  ion  cyclotron 
instability,  the  approximate  altitude  of  this  transition,  and  the 
dependence  of  growth  rate  of  the  ion  cyclotron  instability  on  ambient  high 
latitude  ionospheric  plasma  parameters  is  not  well-defined. 

In  addition  to  observations  indicating  existence  of  EIC  waves  in  the 
bottomside  ionosphere,  the  presence  of  energetic  heavy  ions  at  1400  km 
altitude  was  first  documented  by  Hoffman  et  al.  [1974],  Energy  and  pitch 
angle  distributions  of  various  ion  species  obtained  by  S3- 3  and  PR0GN0Z-7 
satellite  measurements  provide  evidence  for  the  presence  of  energetic 
heavier  ion  species  (N0+  and  0+,  for  example)  at  altitudes  of  a  few 
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thousand  kilometers  or  more  (Hultqvist ,  1983;  Yau  et  al. ,  1984).  The  pitch 
angle  distributions  of  ions  ranged  from  being  field-aligned  (ion  beams)  to 
conical,  in  which  case  the  flux  peaks  at  an  angle  to  the  geomagnetic 
field.  The  pitch  angle  distributions  measured  at  energies  less  than  650  eV 
are  relatively  field-aligned,  while  those  at  1  keV  and  above  have  a  conical 
distribution  with  no  field-aligned  component  [Klumpar  et  al . ,  1984]. 
Recent  data  from  DE-1  also  showed  that  the  high  altitude  Polar  Cap  region 
contains  upflowing  heavy  molecular  ions  such  as  N0+ ,  02*  and  N2+  (Craven  et 
al . ,  1985).  These  data  suggest  that  the  source  for  these  energetic  heavy 
ions  might  be  located  in  the  ionosphere  (Hultqvist ,  1983).  Based  on  our 
linear  theoretical  analysis,  we  conjecture  that  the  current  driven  EIC 
waves  at  the  ion  (heavy)  cyclotron  frequency  could  selectively  accelerate 

4.  +  4. 

the  heavier  ions,  such  as  NO  ,  0 2  and  0  ,  in  the  transverse  direction. 
Thus  the  bottomside  ionosphere  can  act  as  a  source  for  the  heavier  ion 
conics  observed  at  magnetospheric  altitudes. 

In  this  paper  we  study  the  EIC  instability  in  the  bottomside 
collisional  ionosphere.  We  find  that  the  growth  rates  of  the  EIC 
instability  in  the  bottomside  ionosphere  maximize  for  k^p^  »  1  .  In  Sec.  2, 
we  present  the  general  dispersion  relation  for  EIC  waves  in  a  collisional 
plasma  which  we  solve  analytically  both  in  the  collisionless,  weakly 
collisional,  and  strongly  collisional  limits.  We  also  solve  the  full 
dispersion  relation  numerically,  using  the  plasma  density  and  temparature 
profiles  obtained  from  the  Sondrestrom  radar  (Kelly ,  1983).  The  dispersion 
relation  is  solved  for  the  maximally  growing  modes,  for 


1  .4, 


corresponding  to  observed  meter  size  irregularities. 
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In  Sec.  3,  we  compute  expected  nonlinear  saturation  amplitudes  of  the 
electrostatic  ion  cyclotron  instability  in  a  collisional  plasma  based  upon 
Dupree's  resonance  broadening  theory  and  compare  with  experimental 
observations.  Finally  in  Sec.  4,  we  summarize  and  discuss  our  results. 


2.  Linear  Theory 

We  consider  a  homogeneous  plasma  in  a  uniform  magnetic  field  along  the 
z-direction,  3  =*  3  z.  The  plasma  is  collisional  and  the  main  source  of  free 
energy  is  due  to  the  magnetic  field  aligned  currents.  In  this  paper  we 
ignore  inhomogeneous  effects,  e.g.,  density  and  temperature  gradients  and, 
as  a  result,  use  local  theory.  Furthermore,  we  assume  the  scale  height  in 


t- 
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the  high  latitude  ionosphere  to  be  much  larger  than  the  parallel 

wavelengths  involved.  In  our  attempt  to  understand  ion  cyclotron 

instabilities  in  the  bottomside  ionosphere,  we  include  collisional  effects 

and  use  the  general  dispersion  relation  derived  by  Clemmow  and  Dougherty 

(1969).  In  order  to  include  collisional  effects  we  use  a  number  conserving 

BGK  collision  model.  This  model  can  also  be  shown  to  conserve  momentum  if 

the  electrons  are  assumed  to  carry  the  bulk  of  the  current  and  the  ions  are 

considered  to  be  at  rest.  Under  these  conditions,  the  BGK  model  reproduces 

well  electron-neutral  (v  ),  electron-ion  (v  .),  and  ion-neutral  (v.  ) 

en  ei  in 

collisions.  However,  the  BGK  model  does  not  reproduce  well  ion-ion  (v^.) 
or  electron-electron  collisions  and  a  different  collision  model,  e.g., 
Fokker-Planck  approximiation,  must  be  used.  Since  we  will  be  considering 
bottomside  F-region  altitudes  (z  <  300-^00  km),  v..  <  v.  is  satisfied. 
Electron  viscous  effects  (v  )  are  negligible  because  the  electron  Larmor 
radius  is  much  smaller  than  the  perpendicular  wavelengths  of  interest.  In 
section  2.2,  we  do  include  a  phenomenological  model  for  in  our 
numerical  computation  of  growth  rates  and  marginal  stability  criteria  for 
the  electrostatic  ion  cyclotron  instability. 

We  utilize  the  dispersion  relation  derived  by  Clemmow  and  Dougherty 
(1969).  We  list  only  the  salient  features  of  the  derivation  here.  The 
equilibrium  function  is  a  drifting  Maxwellian 


f0j  ■  N0j  <2^j)3/2  *  (\  -  Vdj,2>/Vj21 

where  N„,t  m, ,  T • ,  VH,  and  V*  are  the  density,  mass,  temperature,  drift  and 
uj  J  J  uj  j 

thermal  velocities  of  the  jth  species,  respectively.  The  perturbed 
distribution  function  (where  fj  »  fQj  *  f^)  is  given  as,  using  a  BGK 
collision  model, 
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For  a  low  8  (  oirnkT/B  <<  1 )  plasma  and  perturbations  of  the  form 

f1  «  exp  [ — i ( cut  -  k*x)]  the  dispersion  relation  is  given  as  [Xindel  and 

Kennel ,  1971] 
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where  I  (b)  is  the  modifed  Bessel  function  of  ntn  order  and  Z(£)  is  the 

plasma  dispersion  function,  and  v  =  v,  is  the  collision  frequency  of  the 

J  n 


J  species  with  neutrals  and  qj  is  the  charge  of  the  species  (+  e).  In 
the  case  of  electrons,  Coulomb  collisons  can  also  be  included  and  we 
define  v  »  v  ,  +  v 


The 


effects  of  collisions  are  reflected 
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the  shifted 


frequency  oj  (by  iv^)  and  more  importantly  in  the  additional  term 
proportional  to  in  the  denominator  of  Eq.  (3).  In  fact,  this  term 

gives  rise  to  instability  in  the  collisional  domain.  This  local  dispersion 
relation  is  assumed  to  describe  both  the  collisional  bottomside  and 
collisionless  topside  ionosphere.  We  assume  quasi-neutrality  and  solve 
equation  (3)  analytically  in  the  collisionless,  weakly  collisional,  and 
strongly  collisional  (fluid)  limits. 

2.1  Collisionless  Domain 

By  setting  v.  ■  0  in  equation  (3)  we  immediately  recover  the 
collisionless  dispersion  relation  (Drummond  and  Rosenbluth,  1962;  Kindel 
and  Kennel,  1971).  Furthermore  if  we  set  n  »  0  for  electrons  and  n  =  1  for 
ions  and  expand  the  plasma  dispersion  function  Z(£)  such  that  the  electrons 
are  treated  kinetically  and  the  ions  are  treated  in  the  fluid  limit 
(5g  <<  1  and  >>  1 ,  respectively)  we  obtain  their  results: 


ai^/S^  »  1  +  t  (bi ) 
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where  x  »  T  /T .  and  V, 
e  l  d 


7 .  -  V...  We  note  here  that,  as  will  be  shown  in 
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sec.  2.2,  the  numerical  solutions  of  equation  (3)  agrees  with  equations  (4) 
and  (5)  if  n  =  1  for  ions.  However,  for  n  >  1  a  higher  growth  rate  is 
observed,  which  saturates  for  n  =  4.  Consequently,  the  critical  drift 
velocities  given  in  Kindel  and  Kennel  (1971)  are  overestimated  and  the 
growth  rates  are  underestimated. 

2.2  Weakly  Collisional  Limit 


The  analysis  in  this  limit  was  also  attempted  by  Kindel  and  Kennel 
(1971).  We  expand  the  electron  Z-function  in  the  small  argument 
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Here  we  have  assumed  bg  <<  1  which  gives  (bg)  »  1  .  We  expand  the  ion  Z- 
function  in  the  large  argument  limit,  >>  1.  The  ion  contribution  of 


equation  (3)  is  given  as 
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For  v^/Q^  <  1  the  ion  contribution  becomes 


D.  -  -  f  [u  -  v.2f  )  +  iv.  (1  +  u  f  )] 
l  n  L  in  l  n  J 


where 


f  3  z  - -i! — ; — . 

n  a)  +  nQ  +  iv. 

For  v./(u  <  1  the  dispersion  relation  then  becomes  D  =  D  +  Dj  *  0  or 

1  Cl 

(u  -  k  V  ) 

—  [  1  +  i tt r  - — r. — f-. - £ - y — rf ]  =  "  f  (oj  +  iv.  uv.f  ) 
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From  the  real  part  of  equation,  we  obtain  the  usual  result 


ur  =  ^i^-1  +  T  ri  ^bi  ^ 


and  from  the  imaginary  part  we  obtain 


Y/w  -  -  2  /tt  tI\  (b. 


(co  -  k  V  )  v 

i  - r  .I  ^  -  r  i  +  / tt  t — %- 


u/n.fi  +  -]  (12) 


which  _eads  to  the  collisionless  result  of  Drummond  and  Rosenbluth  (1962) 

and  Kindel  and  Kennel  (1971)  in  equation  (5)  when  vg  =  =  0.  Equation 

(12)  shows  that  the  electron  collisions  have  a  destabilizing  influence 

whereas  the  ion  collisions  are  stabilizing.  In  the  short  wavelength  limit, 

v  /k  7  <<  1  ,  the  ion  collisions  override  the  destabilizing  influence  of 

electron  collisions.  However,  we  show  in  the  next  section  that  the 

electron  collisions  dominate  the  ion  collisions  in  the  fluid 

limit,  v  /k  V  >>  1  . 
e  |  e 

2.3  Strongly  collisional  limit 

In  this  limit  the  electron  mean  freepath  X  <<  X  and  thus  v  /k  V  >> 

e  |  e  |  e 

1.  For  vg/k  Vg  >>  1,  we  can  expand  the  electron  Z-function  in  the  large 
argument  limit  while  still  satisfying  IT  <  <  V  ,  obtaining 

f  (b  ) 

De  =  _l_  - — - - -  (13) 

Ee2  -  p^-2  S  (i  ♦  i/2Ee2)  ro(be) 

*1  e 

The  ion  contribution  is  obtained,  as  before,  by  expanding  the  Z-function  in 
the  large  argument  limit,  ?  >>  1  ,  as 


D  .  a  G  -  T  (b . )  — £- 
i  i  i '  u-n. 


where 


G  =  1  - 


1  *  W 


How,  the  dispersion  relation  becomes 


D  =  D  *  D .  *  0 
e  i 


from  which  the  growth  rate  and  real  frequency  can  be  written  as 


1  *  if  (b. )  (l  -  xG)  (l  - - 


(Q.  -  k  7  ) 

i  id  1 


?  ?  r  fh  1 
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where  a  =  (v^/Q^Kk'Vk^) .  For  b.  <<  1  and  v  / J2  <<  1  equations  (16)  and 

e  e  j_  I  1  e  e 

(17)  reduce  to  that  obtained  by  Chaturvedi  (1976)  who  used  a  fluid 
analysis.  We  see  from  equations  (16)  and  (17)  that  the  real  frequency  is 
not  significantly  different  from  the  collisionless  result  (see  equation  4) 
and  the  growth  rate  is  directly  proportional  to  vg  and  but  has  a  peak 


t-3/2 

V  a  T 
e  e 


for  Coulomb 


in  v  at  v  -  fl  (k  /k  )[r„(b  )]  .  Since  v  «  T  J  for  Coulomb 

e  e  e  |  y  L  0  e  J  ee 

collisions  we  see  from  equation  (17)  that  initially  the  growth  rate 

-7/2 

decreases  with  Tg  as  T  and  for  strong  collisions  the  growth  rate 
7/2  6 

increases  as  indicating  a  minimum  as  a  function  of  t  (for  a  fixed 

TV ) .  This  is  in  sharp  contrast  to  the  collisionless  case  where  the  growth 
rate  increases  as  a  function  of  t  (equation  5). 


2.4  Marginal  stability  analysis 

In  the  collisionless  domain  several  authors  have  estimated  the  drift 
velocity  thresholds  for  the  onset  of  the  ion  cyclotron  instability  (see  for 
example,  Drummond  and  Rosenbluth,  1962;  Kindel  and  Kennel,  1971).  Kindel 


and  Kennel  (1971)  have  examined  briefly  the  effects  of  ion-neutral 
collisions  and  weak  electron-neutral  collisions.  Their  expression  for 
threshold  drift  velocity  shows  that  while  the  electron-neutral  collisions 
could  destabilize  the  EIC  instability,  the  ion-neutral  collisions  strongly 
stabilize  the  EIC  instability. 

In  the  weakly  collisional  limit  v  /k  v  <<  1 ,  we  expand  the  electron 

©  |  6 

Z-function  in  the  small  argument  limit.  The  electron  part  of  the 
dispersion  relation  in  the  above  limit  together  with  the  ion  part  in  the 
large  argument  limit  yield  the  dispersion  relation  given  in  equation 
(10).  The  imaginary  part  of  equation  (10)  yields,  for  Y  =  0, 

f  v,  ( 1  +  O)  f  ) 

u  -  k  V  -  —  (k  V  -  A  V  )  — — -  (18) 

r  '  c  A  •  e  e  (1  -  vff *) 


Where  VQ  is  the  threshold  drift  velocity.  Considering  the  resonance  at  the 
first  harmonic  we  obtain  from  equation  (18) 


V-‘.  - V>.V-V.  'VI-. 
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V  u>  tT,  v/ttv.  u>  r, 

—  >  — —  +  — I  (i  — ® — )(_L)h  +  r  1 )  (19) 

V  k  V  u  k  V  nw-Qn  u-Q  ; 

e  |  e  /ir  j  e 

which  agrees  with  that  of  Kindel  and  Kennel  (1971).  Equation  (19)  shows 
that  in  a  weakly  collisional  limit  the  critical  drift  velocity  is  reduced 
by  electron  collisions  and  increased  by  ion  collisions.  For  a  given  ion 
temperature,  increasing  the  electron  temperature  increases  the  critical 
drift  velocity;  however,  electron  collisions  reduce  the  influence  of  ion 
collisions. 

In  the  strongly  collisional  limit  ve/k^Vg  >>  1  and  >  u>,  the  large 
argument  expansion  of  the  electron  and  ion  Z-functicns  yields  the 
dispersion  relation,  equation  (15),  from  which 


_ W  » ° _ 

r0(t>e)  *  <>)  -  2  UeM.  -  »v0)«V 


"  ♦  r,[-  ‘l)-o 


ID  '  2, 


where  r.tb^  (1  +  tG)  and  u  -  u>  +  iv^.  From  the  imaginary  part  of 

(20),  we  obtain  for  (u  -  k  V  )/k  V  <<  1 

,  c  |  e 


-  —  ■  +  t  r  — *— 

k  V  1  v 

I  e  € 


Vi  ,  “r  /l«il  ! 
- 2 — 2  /ir  ~  e  J 
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where  is  obtained  from  the  real  part  as 


u>r  -  flj  [l  +  Tf/O  ♦  t)1 


Equation  (21  )  shows  that  in  the  strongly  collisional  limit  the  critical 
drift  velocity  has  a  markedly  different  behaviour  than  in  the  weakly 


collisional  or  collisionless  limit. 


Even  for  vi  =  0  the  electron 


collisions  reduce  the  critical  drift  velocity,  which  is  not  the  case  in  the 
weakly  collisional  limit  (equation  19).  Furthermore,  ion  collisions  affect 
the  collisional  result  differently  depending  on  how  strong  the  collision 
frequency  is  compared  to  (w^  -  fl,.),  as  can  be  seen  from  the  second  term 
on  the  right  hand  side  of  equation  (21).  More  importantly,  the  critical 
drift  velocity  has  a  different  behaviour  as  a  function  of  t,  showing  a 
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minimum  in  some  domains,  compared  to  the  collisionless  case,  where  it  was 
shown  to  increase  linearly  with  t  (Drummond  and  Rosenbluth,  1962). 


2.5  Numerical  Results  Applied  to  Auroral  F-region 

The  dispersion  relation  given  in  equation  (3)  describes  the  obliquely 
propagating  ion  cyclotron  instability  and  ion  acoustic  instability  as 
well  (k  -  0,  t  >  1).  In  this  paper,  we  confine  ourselves  to  the  analysis 
of  the  ion  cyclotron  instability.  We  illustrate  the  properties  of  the  EIC 
instability  by  solving  the  dispersion  relation  numerically  using  parameters 


typical  of  the  high  latitude  ionosphere  and  compare  the  numerical  results 

with  the  analytical  results  in  various  regimes  of  collisionality  where 

possible.  We  find  that  the  growth  rate  of  the  collisional  ion  cyclotron 

instability  maximizes  at  k  p.  *  1.  Furthemore,  k  p.  *  1  also  corresponds  to 

1  +  1 

modes  with  X  ■  11m  (pi  -  2.56m  for  0  ),  which  is  on  the  order  of 
typically  observed  scale  size  [Fejer  et  al.,  1984;  Bering,  1984]. 


Typical  bottomside  auroral  zone  ionospheric  electron  density  and 

temperature  profiles  based  on  Sondrestrom  incoherent  scatter  radar  data 

(Kelly,  1983;  C.L,  Rlno,  private  communication)  are  shown  in  Fig.  1  and  2 

respectively.  Curve  B  is  taken  from  Kelly  (1983)  the  data  being  simply 

extrapolated  to  600  km.  The  two  electron  density  profiles  (A  and  C)  are 

expected  to  bracket  both  quiet  and  disturbed  conditions  with  corresponding 

F-peak  densities  of  roughly  10^  and  10^  cm"^,  respectively.  The  electron- 

neutral  (v  )  ,  electron-ion  (v  ,),  ion-ion  (v, .)  and  ion  neutral  (v.  ) 
en  ei  ii  In 

collisions  based  upon  the  data  in  Fig.  1  and  2  are  calculated  using 

standard  expressions  (Spitzer,  1962;  Banks  and  Kockarts,  1973).  Fig.  3 

gives  the  total  electron  collision  frequency  (vg  -  v  +  vgi),  while  Fig.  4 

displays  the  ion  collision  (v.  =  v.  +  v,  . )  frequency.  These  collision 

i  in  ii 

frequencies  are  calculated  using  the  electron  density  profiles  of  Fig.  1 
(Curve  2,  Kelly,  1983).  Neutral  particle  density  profiles  are  taken  from 
Banks  and  Kockarts,  (1973)  and  fall  off  exponentially  with  altitude. 

In  order  to  illustrate  the  destabilizing  effects  of  electron  neutral 
and  electron  ion  collisions  on  the  ion  cyclotron  instability,  we  solve  the 
dispersion  relation  Eq.  (3)  numerically  varying  vg  =  vgn  +  vgiand  keeping 
other  parameters  fixed.  To  illustrate  the  effect  of  collisionality  of  the 
plasma,  Fig.  5  shows  the  growth  rate  of  the  electrostatic  ion  cyclotron 
instability  (Eq.  (3))  as  a  function  of  altitude  for  profile  B  of  Fig.  1, 


assuming  a  pure  0+,  e“  plasma  and  taking  k  /k  -  0.06,  Vd/Ve  -  0.3.  The 

2  2 1  1  u 

growth  rate  shown  is  based  on  b^  =  k^p^  /2  -  1  which  gives  maximum 
growth.  Curve  1  represents  the  collisionless  growth  rate  (v  =  0  =  v^). 

The  decrease  in  growth  rate  is  due  to  the  decrease  in  x  as  one  goes  down  in 


altitude. 


Curve  2  (where  *  0,  =  0),  however,  displays  the  strong 


destabilization  in  the  100  to  150  km  region  and  mild  destabilization  around 


300  km. 


Curve  3,  (where  ve  *  0,  *  0),  on  the  other  hand,  shows  the 


stabilizing  influence  of  ion  collisions.  The  damping  influence  of  ion-ion 

collisions  is  phenomenologically  represented  by  (v^b.^  )/(1  +  bi  ) ,  as 

discussed  in  the  next  paragraph.  As  a  comparison  we  show  in  curve  4  the 

collisionless  growth  with  only  ion  damping  included  (v^  *  0,  v  «  0).  Curve 

3  represents  the  growth  rate  for  completely  realistic  ionospheric 

conditions  for  which  parameters  such  as  plasma  density,  electron  and  ion 

temparatures  and  neutral  density  vary  with  altitude.  Curves  3  and  4 

clearly  establish  the  importance  of  the  electron  collisions,  which  overcome 

the  ion  damping  in  the  region  100  <  z  <  130  and  significantly  enhance  the 

growth  rate  in  the  region  150  <  z  <  300. 

Figure  6  displays  the  growth  rate  of  the  electrostatic  ion  cyclotron 

instability  including  both  electron  and  ion  collisions  for  the  three 

electron  density  profiles  in  Fig.  1  and  temperatures  shown  in  Fig.  2 

assuming  an  0+  ion  plasma.  Although  not  explicitly  calculated  in  Eq.  (3), 

the  ion-ion  collisional  damping  (v^)  is  included  in  the  calculation  of  the 

growth  rates  of  Fig.  6  by  modeling  this  damping  with  the 

expression  -v^k  p ^  / ( 1  +  k  )  (Y.C.  Lee,  private  communication).  This 

expression  reduces  to  the  appropriate  viscous  damping  in  both  the 

magnetized  (k  <  1  )  and  unmagnetized  (k  p^  >>  1)  limits.  Since  we  are 

concerned  with  bottomside  characteristics  of  EIC  waves  (v,,  <  v.  )  ,  the 

11  —  in 

growth  rates  will  not  be  extremely  sensitive  to  v. .  in  any  case.  The  ion- 

11  1/2 

neutral  collisions  were  treated  by  3GK  model  with  v.  <*  n  T.  ,  where  nn 

m  o  l  u 

is  the  neutral  density.  Finally,  for  Fig.  6  the  following  parameters 

k  /k  -  0.06  ,  k  p.  -  /2,  and  VVVQ  =  0.3  are  used.  Curves  A,  B  and  C 
II  1  1  d  e 

correspond  to  the  profiles  given  in  Fig.  1.  It  can  be  seen  from  curve  A 
that  the  growth  rate  maximizes  around  400  km  and  progressively  decreases  at 
lower  altitudes  due  to  increased  damping  by  ion-neutral  and  ion-ion 
collisions  and  due  to  a  decrease  in  the  electron  temperature  (t  decreases 
with  decreasing  altitude;  see  Fig.  2).  Curve  B,  representing  a  typical 


ionosphere  has  smaller  growth  rate  than  curve  A  above  250  km  and  slightly 
higher  growth  rate  than  curve  A  below  250  km.  Although  electron  collisions 
lead  to  higher  growth,  the  decrease  in  t  has  greater  effect  on  the 
collisional  instability  than  on  the  collisionless  instability;  hence,  the 
decrease  in  the  growth  rate  above  250  km.  The  increase  in  the  growth  rate 
at  lower  altitudes  is  mainly  due  to  increased  electron  collisions.  Curve 
C,  representative  of  a  disturbed  ionosphere,  has  peak  growth  at  -  250  km  in 
the  F-region,  the  decrease  in  the  growth  rate  at  higher  altitudes  being 
mainly  due  to  increased  ion-ion  collisional  damping.  Eventually  the  ion 
damping  decreases  at  higher  altitudes  where  the  collisionless  theory  is 
valid  above  which  altitude  the  growth  rate  increases  due  to  increasing  t. 


Figure  7  gives  the  threshold  drift  velocities  (Vd/Ve)  for  the  EIC 
instability  for  the  three  electron  density  profiles  of  Fig.  1  assuming  an 
0+  ion  plasma  taking  k^/k^  ”  and  for  maximally  growing  modes 
with  k  =*  v^2  .  For  profile  B,  the  electrostatic  ion  cyclotron 
instability  could  be  destabilized  in  the  bottomside  ionosphere  (  <  300  km) 
for  current  drift  velocities  on  the  order  of  a  tenth  of  electron  thermal 


velocites,  Vd  i  0.1  Vg.  Curves  A,B  and  C  again  are  obtained  using  the 
density  profiles  given  in  Fig.  1.  Comparing  thi3  figure  with  Fig.  6,  we 
see  that  lower  drift  velocities  are  '*«•:' -id  at  larger  growth  rates.  This 
figure  also  shows  that  the  collisional  EIC  instability  can  al3o  be  excited 
in  the  E  region,  although  it  shows  that  it  is  easier  to  excite  EIC  in 
moderately  collisional  domain  (Curve  A)  than  under  disturbed  conditions 
(Curve  C). 


Figure  8  shows  the  growth  rate  of  the  EIC  instability  including  both 


electron  and  ion  collision  for  the  three  electron  density  profiles  in  Fig  1 
assuming  an  N0+  ion  plasma.  The  growth  rates  of  the  EIC  instability  for  an 
N0+  plasma  are  similar  to  (though  slightly  lower  than)  the  growth  rate  for 
the  0+  plasma  (Fig.  6).  This  is  consistent  with  Eqs.  (5),  (12)  and  (17) 


which  indicate  that  the  Y/Q^  is  not  a  strong  function  of  ion  mass.  Note 
that  the  growth  rate  is  normalized  to  respective  gyro  frequencies  in  each 


case. 


Fig.  9  displays  the  threshold  drift  velocities  (Vc/Ve)  for  the  EIC 
instability  including  both  electron  and  ion  collisions  using  the  density 
profiles  of  Fig.  1  assuming  N0+  ion  plasma.  Comparing  Fig.  9  with  Fig.  7 
(0+  plasma)  one  notes  a  slight  decrease  in  the  threshold  drift  velocity  as 
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one  goes  from  0+  to  N0+.  This  is  also  consistent  with  previous  studies 
(Drummond  and  Rosenbluth,  1962,  Kindel  and  Kennel,  1971)  which  have  shown 
that  increasing  the  ion  mass  decreases  the  threshold  current  drift  velocity 
at  all  Tg/Ti  to  excite  the  EIC  instability.  We  have  also  calculated  the 
growth  rates  and  the  critical  drift  velocities  for  the  Hydrogen  plasma.  We 
find  the  growth  rates  of  H+  EIC  waves  to  be  smaller  and  the  critical  drift 
velocities  for  excitation  of  H+  EIC  waves  to  be  much  higher;  for  instance, 
Vc/7 e  »  0.4  for  H+  EIC  wave  excitation  for  the  typical  parameters  used  for 
0 *  calculations,  as  compared  to  Vc/V  -  o.l  for  0+  excitation.  These 
critical  drift  velocities  correspond  to  extremeley  large  currents  and  thus 
we  can  conclude  that  the  H+  EIC  waves  are  not.  likely  to  be  excited  in  the 
bottomside  ionosphere.  For  this  reason  we  do  not  present  the  detailed 
results  of  H+  EIC  wave  calculations. 


To  further  emphasize  the  heavy  ion  EIC  excitation  in  the  F  region,  we 
show  the  threshold  currents  in  the  altitude  region  100  km  to  600  km.  We 
use  the  density  profile  (ng)  from  Curve  A  of  Fig.  1,  Te  from  Fig.  2  and 
( Vc/Ve)  from  Figs.  7  and  9  to  calculate  the  threshold  current  J  -  nee'I  . 
Figure  10  shows  the  threshold  currents  required  to  excite  0+  EIC  and  N0+ 
EIC  in  0+  and  N0+  plasmas,  respectively  for  k  -  /2  and  k^/k  -  0.06.  It 
is  clear  from  Fig.  10  that  N0+  EIC  has  lower  thresholds  than  0+  EIC  and 
the  altitude  region  between  150  km  and  200  km  and  the  region  above  400  km 
have  lower  thresholds.  The  currents  shown  near  100  km  may  not  be  accurate 
because  the  theory  may  not  be  valid  in  this  region  due  to  demagnetization 


of  the  ions  (v.  /n.  -  1 ) 
in  i 


The  critical  currents  corresponding  to  curves  3 


of  Figs.  1,  7  and  9  (not  shown  here)  are  of  the  order  200  uA/m  -  1000 
2 

uA/m  .  Thus,  for  typical  ionospheres  such  as  those  given  by  curve  B,  one 
requires  large  currents  to  support  0+  or  N0+  EIC.  Although  our 
calculations  are  applicable  to  plasma  with  single  ion  species,  based  on  the 
results  presented  in  Figs.  7,  9,  and  10  we  can  draw  qualitative  conclusions 
regarding  a  multispecies  plasma,  namely  that  when  H+,  0+ ,  0*  and  N0+ 
species  are  present,  EIC  waves  of  heavier  species  are  more  susceptible  to 
destabilization  by  parallel  currents.  Thus  we  conclude  that  N0+  or  0*  EIC 
instability  can  be  excited  more  easily  in  the  bottomside  F  region  than  the 
H  +  or  0+  EIC  instability.  This  is  consistent  with  the  photometric 
observations  of  Martelli  et  al.  [1971]. 


3.  Nonlinear  Saturated  Amplitudes 

In  the  previous  section  we  showed  that  the  collisional  electrostatic 
ion  cyclotron  instability  could  be  important  in  the  bottomside  high 
latitude  ionosphere.  In  this  section  we  discuss  a  possible  nonlinear 
mechanism  that  might  saturate  these  modes.  In  the  literature  several 
nonlinear  saturation  mechanisms  for  the  electrostatic  ion  cyclotron 
instability  have  been  studied.  Davidson  (1972)  has  discussed  large 

amplitude  trapping  effects  which  become  important  when  u  x  >>  1  where 

1/2  ^  G 
a).  =  k  (e$/m  )  is  the  electron  trapping  frequency  and  x  =  k  Av  is  the 
tie  c  | 

wave  correlation  time,  <j>  is  the  wavepotential ,  _k  is  the  wavevector, 
and  AVp  is  the  spread  in  wave  phase  velocities  parallel  to  the  magnetic 
field.  Collisional  effects  will  prevent  strong  electron  trapping.  Another 
mechanism  is  ion  resonance  broadening  [Dum  and  Dupree,  1970].  In  this  work 
we  extend  Dum  and  Dupree’s  theory  to  the  case  of  the  collisional  EIC 
instability.  Physically,  the  effect  of  the  turbulence  generated  by  the 
instability  may  be  interpreted  as  an  enhanced  ion  viscosity.  In  this 
section  we  shall  present  only  estimates  of  the  saturation  amplitudes  for 
typical  ionospheric  parameters  and  not  a  comprehensive  full  F-region  study. 

The  nonlinear  dispersion  relation  for  the  electrostatic  ion-cyclotron 
modes,  including  the  effects  of  resonance  broadening,  may  be  obtained  by 
substituting  w  a  cu  +  iAa)W  in  lieu  of  oj  in  the  ion  contribution  to  Eq.  (3)  [ 
Dum  and  Dupree,  1970],  Here  A<uW  is  the  modification  introduced  by  the 
effect  of  a  spectrum  of  finite  amplitude  modes  on  the  particle  orbits 
(resonance  broadening  effect).  Thus, 


1  +  e^k,  <*>)  +  a)  +  iA  ujW) 


is  the  nonlinear  dispersion  relation  for  the  EIC  modes,  where  e.  is  the 

n  q  ^ 

linear  dielectric  constant  and  ej  Is  the  nonlinear  dielectric  constant. 
We  note  that  Dum  and  Dupree  [1970]  have  shown  the  effect  of  resonance 
broadening  is  important  only  for  ions  in  the  EIC  wave  case.  Hence  the 
contribution  coming  from  the  resonance  broadening  effect  on  electrons  has 
been  ignored  in  the  expression  (23).  One  may  obtain  an  estimate  for  the 
nonlinear  saturated  amplitudes  of  EIC  modes  from  (23)  by  setting  the 
nonlinear  growth  rate  Y 
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Wl  A  w 

y,  -  Aw 

k 


where  is  the  linear  growth  rate  of  the  mode  with  wavenumber  k.  For  the 

w 

EIC  instability  case,  Aw  may  be  approximated  as  [Dum  and  Dupree,  1970], 


/  •  U-  0JK-5)2  - 1]' 


Here  n^  is  the  saturated  nonlinear  amplitude  of  the  wave,  and,  n^  is  the 
threshold  level  of  the  nonlinear  amplitude  at  which  the  resonance 
broadening  effects  assume  importance.  The  expression  for  n1?  is 


n,  T.  0.  „ 

J_  =  ( — L)( — i_)2 

n  lZT  Mk  v.J 


( w,  "  0 .  ) 
v  k  1  ‘ 


e  j_  i  aitF1  Cp± )  1 


T  w  -  Q 

[—1]  - 3 -  f - i) 

ZTe  u^F^172  ^  fli 


where 


F^x) 


{b  -  i  (§)  +  5  (f )  ♦  for  x  S  2 

~  x  2  3 


It  may  be  noted  that  the  broadening  of  resonances  above  corresponds  to 
an  enhancement  in  the  group  of  ions  that  can  exchange  energy  with  the  wave 
via  resonance  interaction.  Stabilization  of  the  mode  results  when  the 
resonance  is  broadened  to  an  extent  such  that  the  bulk  of  the  ion 
distribution  interacts  with  the  wave  and  absorbs  energy  from  it,  thereby, 
leading  to  a  steady  finite  amplitude  nonlinear  state.  We  now  use 
expressions  (24)-(25),  to  compute  the  nonlinear  saturated  amplitudes  for 
the  collisions!  EIC  modes.  From  (24)  and  (25),  one  has 


(u-  OjHt-i)2  -  .] 


which  leads  to 


K)2 


(— )2  ■  r— )2r  1  *  — s — 1 

'•n  ’  1  n J  L  ,  „  a2j 


16 


or 


(-L.  2  _ i_  4  (JL_ i  2 

ni  o  ?  U  V.  1  {  Q.  1 

f_L)2  -  _® _ I  L  1 _ 

\n  ;  k  v 

F  (— 1 — L) 

I1-  Q  J 
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[i 


— —  f  k  V 
,  2  2  1  id 
k  v  * 

1  e  - 


u  )  -  k^pfv.  -v. 
r  ri  xi  in' 


UP  " 


(29) 


using  (17)  and  including  collisions, 


For  the  case,  Te  -  Ti ,  one  has 


U  -  n.)  -  n.r, 


which  leads  to 


Q-Y 


lk  v. ;  1 

1  1  [1 


k  v 


M  Q 


[-T2  ^d 

k  v  * 

I  e 


,  2  2 

k  p . v. . - 
1  I  i  ii 
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}] 


(30) 


The  expression  (30)  gives  the  nonlinear  saturated  amplitudes  of  density 
fluctuations  associated  with  the  collisional  ion-cyclotron  modes  due  to 
resonance  broadening  effects  in  the  presence  of  both  electron  Ue)»  ion- 
ion  (vu)  and  ion-neutral (\k  )  collisions  applicable  to  bottomside 
ionospheric  altitudes.  A  comparison  with  the  collisionless  case,  discussed 
by  Dum  and  Dupree  Cl 970],  shows  that  the  saturated  amplitudes  can  be  higher 
in  the  collisional  case,  by  an  increment  of  approximately  -  tr[(v./V  ) 

(l/k  \  )]  .  We  note  that,  in  the  above  comparison,  it  is  implied  that 

I  e 

Ykc/(^k  "  <<  1 »  where  is  the  collisionless  growth  rate. 


U.  Summary  and  Discussion 

In  this  study  we  have  presented  a  comprehensive  treatment  of  the 
theory  of  the  electrostatic  ion  cyclotron  instability  in  the  collisional 
bottomside  ionosphere.  We  have  derived  the  linear  dispersion  relations  for 
the  ion  cyclotron  instability  including  both  electron  and  ion  collisions 
and  solved  analytically  the  dispersion  relation  in  three  limits: 
collisionless,  weakly  collisional,  and  collisional.  In  addition,  we  have 
solved  the  linear  dispersion  relation  exactly  using  numerical  methods  using 
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parameters  typical  of  the  high  latitude  bottomside  ionosphere.  We  note 
that  the  linear  growth  rates  in  the  bottomside  collisional  regime  can  be  of 
the  same  order  as  the  corresponding  topside  collisionless  growth  rates.  We 
have  considered  three  cases  corresponding  to  three  different  ionospheric 
conditions.  The  normal  ionosphere  is  bracketed  by  two  others  with  minimum 
density  around  1014  cnT^  at  the  lower  end  and  maximum  peak  density  of  10^ 

cm"  3  at  the  high  end.  The  former  represents  a  weakly  collisional 
-4 

case  (v  /ft  <10  )  and  this  can  be  considered  as  an  extension  of  the 
e  e 

previous  collisionless  EIC  instability  theories,  the  latter  represents 
disturbed  conditions,  yielding  high  electron  collisions  and  ion 
collisions.  We  have  also  considered  pure  0+  and  pure  M0+  plasmas  and 
solved  the  dispersion  relation  for  single  species  in  both  cases.  We  arrive 
at  the  following  conclusions  based  on  the  linear  calculation. 


a.  Given  the  realistic  ionospheric  parameters,  we  find  that  the  theory  of 
collisionless  EIC  instability  (Kindel  and  Kennel,  1971 )  can  not  be 
extended  to  the  bottomside  F  region  (z  <  400  km).  Since,  the  critical 
drift  velocities  and  growth  rates  are  sensitve  to  n0  and  Tg/T^,  a 
theory  of  collisional  EIC  using  self-consistent  collision  frequencies 
and  Ta/T ^  should  be  used  for  accurate  prediction  of  EIC  excitation. 


b.  Electron  collisions  destabilize  the  EIC  wave  which  results  in  higher 
growth  rates.  However,  the  ion-ion  collisions  damp  the  modes,  the 
damping  being  equal 

The 


the  k^p^  <<  1  limit. 

for  k  p,  >  1  . 

1  A  ~ 


to  v, ,  in  k  p .  >>  1  limit 
ii  /i 

EIC  instability  growth  rate  maximizes 


2  2 

and  v. . k  p .  /2  in 
ii 


c.  For  realistic  profiles  of  the  ionosphere,  the  collisional  EIC  wave 
growth  rates  are  reduced  by  ion-neutral  collisions  in  the  bottomside  F 
region  and  by  ion-ion  collisions  above  300km,  as  expected. 


d.  For  profiles  given  in  Fig.  1,  the  collisionless  theory  (ion  collisions 
dominate  electron  collisions)  prevails  for  curve  A  and  collisional 
theory  (electron  collisions  are  stronger  and  dominate  ion  collisions) 
prevails  for  curve  C. 


e.  The  growth  rate  maximizes  in  the  weakly  collisional  (curve  A,  Fig.  1) 
regime  at  an  altitude  of  400  km  and  in  the  strongly  collisional  (curve 
C,  Fig.  1)  regime  at  an  altitude  of  250  km.  In  the  collisional 
domain,  (curve  C,  Fig  6)  the  growth  rate  has  minimum  at  an  altitude  of 
300  km  and  achieves  a  further  lower  value  at  E  region  altitudes  (150 
km  or  less). 

f.  The  critical  drift  velocity  required  to  excite  the  EIC  waves  is 

sensitive  to  the  electron  densites  and  thus  are  quite  different  in  the 
collisionless  and  collisional  domains.  It  is  smaller  at  higher 

altitudes,  z  <  400  km,  and  larger  at  lower  altitudes,  z  <  200  km. 

g.  For  realistic  ionospheric  profiles,  the  growth  rate  is  weakly 

dependent  on  the  mass  ratio  (m0/m^).  However,  we  find  that  the 

critical  drift  velocity  is  much  smaller  for  EIC  waves  corresponding  to 
gyrofrequencies  of  heavier  ions  than  the  critical  drift  corresponding 
to  lighter  ion  (H+)  EIC  wave  excitation.  Our  single  species 
calculation  shows  that  N0+  or  02+  EIC  instability  are  much  easier  to 
excite  than  0+  or  H+  EIC  instability.  Thus,  in  a  multi-component 
plasma  EIC  waves  at  N0+  will  be  first  excited  before,  for  example,  0+ 
EIC  waves  and  in  the  bottomside  ionosphere  EIC  instability  at  H+ 

gyrofrequency  is  unlikely  to  be  excited. 

In  addition,  we  have  discussed  a  nonlinear  saturation  mechanism  for 
the  electrostatic  ion  cyclotron  instability  in  the  collisional  regime  based 
upon  Dupree's  ion  resonance  broadening  theory.  Dupree's  resonance 
broadening  theory  is  valid  only  for  a  broad  spectrum  of  waves  which  are 
assumed  isotropic  and  incoherent,  i.e.,  the  wave  coherence  time  (spectral 
width)  is  much  shorter  than  an  ion  cyclotron  period.  For  a  narrow 
(coherent)  wave  spectrum,  resonance  broadening  theory  must  be  replaced  by 
an  analysis  which  considers  exact  ion  dynamics  in  a  coherent  wave 
[Davidson,  1972;  Aamodt,  1970].  Recent  nonlinear  numerical  simulations 
[Pritchett  et  al.,  1981;  Okuda  and  Ashour-Abdalla,  1981]  of  the 
collisionless  ion  cyclotron  instability  indicate  a  large  amplitude  coherent 
wave  spectrum  in  the  nonlinear  regime.  However  these  simulations  also  show 
several  modes  excited  in  the  saturation  phase  of  the  ion  cyclotron 


instability.  To  our  knowledge,  there  have  been  no  nonlinear  simulations  of 
the  collisional  current-driven  ion  cyclotron  instability.  We  shall  now 
discuss  some  particular  examples  of  the  auroral  ionospheric  application  of 
the  collisional  EIC  instability. 


In  a  recent  study  using  sounding  rockets,  Yau  et  al.  (1983)  have 
recorded  0*  waves  and  strong  perpendicular  ion  acceleration  (ion  conics) 
events  at  low  altitudes  (400-600  km)  in  the  high  latitude  ionosphere  during 
the  expansive  phase  of  an  auroral  substorm.  Large  relative  density 
fluctuations  (<  10$)  were  also  recorded  near  the  ion  conic  acceleration 
regions.  One  major  candidate  far  conic  formation  is  EIC  wave  heating 
[Ungstrup  et  al . ,  1979;  Okuda  and  Ashour-Abdaila,  1981].  Considering 


typical  high  latitude  parameters,  i.e.,  T  *  T. 
corresponding  to  curve  A  of  Fig.  1,  (n  -  3  x  10 


0.1  eV,  electron  density 

v  /£2  -  2  x  10"5, 

e  e 


v./Q.  -  7  x  10  the  growth  rate  of  0*  EIC  waves,  with  X  -  11m  fk  p.  - 

li  i  v  i  i 


1.4,  -  2.56m  for 


1 

0  ions)  and  X  -  200  m,  is  0.07  ft.  or  21 
+  I  1 

taking  ft^  for  0  in  a  0.5  Gauss  magnetic  field  as  300  Hz  (Fig.  6,  curve 

A).  The  corresponding  growth  rate  for  NO*  EIC  waves  (X  -  15m)  is 

0.08  Sh  or  12  sec'1,  taking  ft.  to  be  150  Hz.  The  critical  drift  velocities 

can  be  read  off  Fig.  7  for  0*  ions.  For  the  above  parameters  the  critical 

drift  velocity  at  an  altitude  of  300  km,  for  example,  is  0.09  V  (Curve  A, 

For 


sec 


quiet  conditions)  which  corresponds  to  53  pA/m2  for  n  -  3  x  1010/m^. 

NO*  ions,  under  the  same  conditions,  the  critical  drift  velocity 


is  -  0.06  V  (Fig.  9)  which  corresponds  to  40pA/m  .  For  other  profiles, 

Curves  B  (normal  conditions)  and  C  (disturbed  conditions)  of  Fig  1  critical 

2 

drift  velocities  are  higher  and  require  currents  >  100  pA/m  and 
2 

-  lOOOpA/m  ,  respectively.  In  view  of  these  results  we  find  that  the 
ionospheric  conditions  corresponding  to  curve  A  of  Fig.  1  are  most 
conducive  to  EIC  excitation  at  NO  and  0^  gyrofrequencies .  From  equation 
(30),  the  saturated  amplitude  for  typical  parameters  i.e., 


T. 


1 0_il  sec"1 


■0.1  eV,  ft.  *  3  x  102  sec"1,  v 

-3  1  -2  6 

10  ,  v-^/ft^  -  10  is  found  to  be  approximately  n.j/n0  *  0.24 


v.  /ft. 
in  l 

according  to  ion  resonance-broadening  theory  for  k  p  »  1.2,  k  /k  =»  0.1, 

6  1  1  II 

and  V  >  2  x  10  cm/sec.  We  note  that  the  saturation  amplitude  n^nQ  of 

EIC  waves  in  the  collisionless  regime,  using  Eq.  (?0),  is  slightly 


smaller 


n  /n 
-10 


.22)  than  the  collisional  case. 


Compelling  evidence  of  the  presence  of  energetic  heavy  ionspheric  ions 

at  high  altitudes  (>  1000  km)  (Hultqvist,  1983;  Klumpar  et  al . ,  1 98 M ;  Yau 

et  al . ,  1984)  leads  one  to  suggest  that  the  bottomside  ionosphere,  which  is 

a  source  of  heavy  ions,  is  indeed  acting  as  an  accelerating  region  for 

these  heavy  ions  (Hultqvist ,  1983).  Furthermore,  energetic  heavy  ions,  N0+ 
+ 

and  0  ,  have  been  observed  at  outer  magnetospheric  altitudes  on  days  of 
high  geomagnetic  activity  (Klecker  et  al.,  1985;  Craven  et  al . ,  1985). 
High  geomagnetic  activity  at  higher  altitudes  is  usually  indicative  of  a 
disturbed  ionosphere,  which  could  sustain  large  ionospheric  currents  that 
drive  collisional  SIC  waves  in  the  bottomside  ionospher  unstable.  Thus,  we 
conjecture  that  the  ion  cyclotron  instability  at  the  ion  (heavy)  cyclotron 
frequency  ,  driven  unstable  by  parallel  field-aligned  currents  may 
selectively  accelerate  the  N0  +  ,  0^  and  0+  ions  in  the  transverse  direction 
at  low  altitudes  and  these  transversely  accelerated  ions  form  the  heavy  ion 
conics  at  the  higher  altitudes. 
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Typical  ionospheric  electron  density  profiles  based  on  Sondrestrom 
incoherent  scatter  radar  data.  Curve  A  represents  a  quiet  condition, 
curve  B  an  intermediate  electron  density,  and  curve  C  a  disturbed 
condition. 
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Fig.  5  Growth  rate  Y/S^  of  electrostatic  ion  cyclotron  instability  vs. 

altitude  for  0+  ion  plasma  using  curve  B  of  Fig.  1  and  taking  i  as 

varying  with  altitude  and  vg  *  ”  0  (curve  1);  vg  -  vg(z),  vi  =*  0 

(curve  2);  v  -  v  (z),  v.  •  v.(z)  (curve  3);  and  v  *■  0,  v.  -  v.(z) 

6  6  11  6  11 

(curve  iJ).  Here  k  /k  «  0.06,  k  p.  -  Jz  and  VH/Va  -  0.3. 
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electrostatic  ion  cyclotron  instability  vs.  altitude  for  curves  A-C  of 
Fig.  1  for  0+  ion  plasma  and  k  /k  -  0.06,  k  p.  -  fl. 
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jctrostatic  ion  cyclotron  instability  vs. 
asma  using  curves  A-C  of  Fig.  1  with  fixed 
k  p.  -  /2  and  VH/V  -  0.3. 
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Fig.  9  Critical  current  drift  velocity  (Vc/’/e)  ^or  excitation  of 

electrostatic  ion  cyclotron  instability  vs.  altitude  for  curves  A-C  of 

Fig.  1  for  N0+  ion  plasma  and  k  /k  -  0.06,  k  p  »  /2. 
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J  =  neeVd  (/xAmp/m2) 

10  Threshold  currents  in  yA/m c  required  to  excite  0+  and  N0+  EIC  waves  in 
0+  and  N0+  plasmas,  respectively.  The  critical  currents  are 
calculated  for  J  -  ng  e  Vd  using  Vc/Ve  from  curves  A  of  Figs.  7  and  9, 
ng  from  curve  A  of  Fig.  1  and  Tg  from  Fig.  2.  The  parameters  used  are 
ko,  -  /2  and  k  /k  -  0.06. 
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SYRACUSE,  NY  13201 

01 CY  ATTN  G.  MILLMAN 

GEOPHYSICAL  INSTITUTE 
UNIVERSITY  QF  ALASKA 
FAIRBANKS,  AK  99701 

(ALL  CLASS  ATTN:  SECURITY  OFFICER) 
01 CY  ATTN  T.N.  DAVIS  (UNCLASS  ONLY) 
0 1 CY  ATTN  TECHNICAL  LIBRARY 
01 CY  ATTN  NEAL  BROWN  (UNCLASS  ONLY) 

GTE  SYLVANIA,  INC. 

ELECTRONICS  SYSTEMS  GRP-EASTERN  DIV 
77  A  STREET 
NEEDHAM,  MA  02194 

0 1 CY  ATTN  DICK  STEINHOF 


HSS,  INC . 

2  ALFRED  CIRCLE 
BEDFORD,  MA  01730 

0 1 C  Y  ATTN  DONALD  HANSEN 

ILLINOIS,  UNIVERSITY  OF 
107  COBLE  HALL 
150  DAVENPORT  HOUSE 
CHAMPAIGN,  IL  61820 

(ALL  CORRES  ATTN  DAN  MCCLELLAND) 
0 1 C Y  ATTN  K.  YEH 

INSTITUTE  FOR  DEFENSE  ANALYSES 
1801  NO.  BEAUREGARD  STREET 
ALEXANDRIA,  VA  22311 
01 CY  ATTN  J.M.  AEIN 
01 CY  ATTN  ERNEST  BAUER 
01 CY  ATTN  HANS  WOLFARD 
0 1 CY  ATTN  JOEL  BENGSTON 

INTL  TEL  &  TELEGRAPH  CORPORATION 
500  WASHINGTON  AVENUE 
NUTLEY,  NJ  07110 

01 CY  ATTN  TECHNICAL  LIBRARY 

JAYCOR 

11011  TORREYANA  ROAD 

P.O.  BOX  85154 

SAN  DIEGO,  CA  92138 

0 1 CY  ATTN  J.L.  SPERLING 

JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  MD  20810 

0 1 CY  ATTN  DOCUMENT  LIBRARIAN 
01 CY  ATTN  THOMAS  POTEMRA 
01 CY  ATTN  JOHN  DASSOULAS 

KAMAN  SCIENCES  CORP 
P.O.  BOX  7463 

COLORADO  SPRINGS,  CO  80933 
0 1 CY  ATTN  T.  MEAGHER 

KAMAN  TEMPO-CENTER  FOR  ADVANCED 
STUDIES 

816  STATE  STREET  (P.O  DRAWER  QQ) 
SANTA  BARBARA,  CA  93102 
0 1 C Y  ATTN  DAS  I  AC 
0 1 C Y  ATTN  WARREN  S.  KNAPP 
0 1 C Y  ATTN  WILLIAM  MCNAMARA 
0 1 C Y  ATTN  B.  GAMBILL 


LINKAB IT  CORP 
10453  ROSELLE 
SAN  DIEGO,  CA  92121 

01 CY  ATTN  IRWIN  JACOBS 

LOCKHEED  MISSILES  4  SPACE  CO.,  INC 
P.O.  BOX  504 
SUNNYVALE,  CA  94088 
01  CY  ATTN  DEPT  60-12 
01 CY  ATTN  D.R.  CHURCHILL 

LOCKHEED  MISSILES  &  SPACE  CO.,  INC. 
3251  HANOVER  STREET 
PALO  ALTO,  CA  94304 

01 CY  ATTN  MARTIN  WALT  DEPT  52-12 
01 CY  ATTN  W.L.  IMHOF  DEPT  52-12 
01 CY  ATTN  RICHARD  G.  JOHNSON 
DEPT  52-12 

01 CY  ATTN  J.B.  CLADIS  DEPT  52-12 

MARTIN  MARIETTA  CORP 
ORLANDO  DIVISION 
P.O.  BOX  5837 
ORLANDO,  FL  32805 

01 CY  ATTN  R.  HEFFNER 

M.I.T.  LINCOLN  LABORATORY 
P.O.  BOX  73 
LEXINGTON,  MA  02173 

01 CY  ATTN  DAVID  M.  TOWLE 
0 1 C Y  ATTN  L.  LOUCHLIN 
0 1 CY  ATTN  D.  CLARK 

MCDONNEL  DOUGLAS  CORPORATION 
5301  BOLSA  AVENUE 
HUNTINGTON  BEACH,  CA  92647 
01 C  Y  ATTN  N.  HARRIS 
0 1 C Y  ATTN  J.  MOULE 
01 CY  ATTN  GEORGE  MROZ 
01 CY  ATTN  W.  OLSON 
01 CY  ATTN  R.W.  HALPRIN 
0 1  CY  ATTN  TECHNICAL 

LIBRARY  SERVICES 


MISSION  RESEARCH  CORPORATION 
735  STATE  STREET 
SANTA  BARBARA.  CA  93101 
0 1 C  Y  ATTN  P.  FISCHER 
0 1 C Y  ATTN  W.F.  CREVIER 
01 CY  ATTN  STEVEN  L.  GUTSCHE 
0 1 C Y  ATTN  R.  BOGUSCH 
01 CY  ATTN  R.  HENDRICK 
0 1 C Y  ATTN  RALPH  KILB 
01 CY  ATTN  DAVE  SOWLE 
0 1 CY  ATTN  F.  FAJEN 
01CY  ATTN  M.  SCHEIBE 
0 1 CY  ATTN  CONRAD  L.  LONGMIRE 
01 CY  ATTN  B.  WHITE 
0 1 CY  ATTN  R.  STAGAT 

MISSION  RESEARCH  CORP. 

1720  RANDOLPH  ROAD.  S.E. 
ALBUQUERQUE,  NEW  MEXICO  87106 
01 CY  R.  STELLINGWERF 
01 CY  M.  ALME 
0 1 CY  L.  WRIGHT 

MITRE  CORP 

WESTGATE  RESEARCH  PARK 
1820  DOLLY  MADISON  BLVD 
MCLEAN,  VA  22101 
0 1 CY  ATTN  W.  HALL 
0 1 C Y  ATTN  W.  FOSTER 

PACIFIC-SIERRA  RESEARCH  CORP 
123^0  SANTA  MONICA  BLVD. 

LOS  ANGELES,  CA  90025 

0 1 C Y  ATTN  E.C.  FIELD,  JR. 


PHOTQMETRICS ,  INC. 

4  ARROW  DRIVE 
WOBURN,  MA  01801 

0 1 CY  ATTN  IRVING  L.  KOFSKY 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  3027 
BELLEVUE,  WA  98009 

0 1 CY  ATTN  E.J.  FREMOUW 


R  4  D  ASSOCIATES 
P.O.  BOX  9695 
MARINA  DEL  REY,  CA  90291 
0 1 CY  ATTN  FORREST  GILMORE 
01 CY  ATTN  WILLIAM  B.  WRIGHT, 
01 CY  ATTN  WILLIAM  J.  KARZAS 
01 CY  ATTN  H.  ORY 
01 CY  ATTN  C.  MACDONALD 

RAND  CORPORATION,  THE 
1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 
01 CY  ATTN  CULLEN  CRAIN 
01 CY  ATTN  ED  BEDROZI AN 

RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUDBURY,  MA  01776 

01  CY  ATTN  BARBARA  ADAMS 

RIVERSIDE  RESEARCH  INSTITUTE 
330  WEST  42nd  STREET 
NEW  YORK,  NY  10036 

01 CY  ATTN  VINCE  TRAPANI 

SCIENCE  APPLICATIONS,  INC. 

1150  PROSPECT  PLAZA 
LA  JOLLA,  CA  92037 

01 CY  ATTN  LEWIS  M.  LINSON 
0 1 CY  ATTN  DANIEL  A.  HAMLIN 
01  CY  ATTN  E.  FRIEMAN 
0 1 C Y  ATTN  E.A.  STRAKER 
01 CY  ATTN  CURTIS  A.  SMITH 

SCIENCE  APPLICATIONS,  INC 
1710  GOODRIDGE  DR. 

MCLEAN,  VA  22102 
0 1 C Y  J.  COCKAYNE 
01CY  E.  HYMAN 


PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
318  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 
(NO  CLASS  TO  THIS  ADDRESS) 

0 1 C Y  ATTN  IONOSPHERIC  RESEARCH  LAB 


PHYSICAL  DYNAMICS,  INC. 
P.O.  BOX  10367 
OAKLAND,  CA  94610 
ATTN  A.  THOMSON 


SRI  INTERNATIONAL 
333  RAVENSWOOD  AVENUE 
MENLO  PARK,  CA  94025 


0 1  CY 

ATTN 

J.  CASPER 

01  CY 

ATTN 

DONALD  NEILSON 

0 1 C  Y 

ATTN 

ALAN  BURNS 

0  1  CY 

ATTN 

G.  SMITH 

0 1  CY 

ATTN 

R.  TSUNODA 

0 1  C  Y 

ATTN 

DAVID  A.  JOHNSON 

01  CY 

ATTN 

WALTER  G.  CHESNUT 

0 1  CY 

ATTN 

CHARLES  L.  RINO 

0 1 C  Y 

ATTN 

WALTER  JAYE 

01  CY 

ATTN 

J.  VICKREY 

0 1 C  Y 

ATTN 

RAY  L.  LEADABRAND 

0 1  C  Y 

ATTN 

G.  CARPENTER 

01  CY 

ATTN 

G.  PRICE 

0 1  C  Y 

ATTN 

R.  LIVINGSTON 

0 1  CY 

ATTN 

V.  GONZALES 

0 1  C  Y 

ATTN 

D.  MCDANIEL 

TECHNOLOGY  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
BEDFORD,  MA  01730 

0 1 CY  ATTN  W.P.  BOQU 1ST 

TOYON  RESEARCH  CO. 

P.O.  Box  6890 

SANTA  BARBARA,  CA  93111 


0  1  C  Y 

ATTN  JOHN  ISE, ' JR. 

0 1  CY 

ATTN  JOEL  GARBARINO 

TRW  DEFENSE  &  SPACE  SYS  GROUP 

ONE  SPAC 

E  PARK 

REDONDO 

3EACH ,  CA  90278 

0 1  CY 

ATTN  R.  K.  PLEBUCH 

0 1  CY 

ATTN  S.  ALTSCHULER 

0 1  CY 

ATTN  D.  DEE 

01  CY 

ATTN  D/  STOCKWELL 

S NTF/1 575 

VISIDYNE 

SOUTH  BEDFORD  STREET 
BURLINGTON,  MASS  01803 
01 CY  ATTN  W.  REID Y 
0 1 CY  ATTN  J.  CARPENTER 
01CY  ATTN  C.  HUMPHREY 

UNIVERSITY  OF  PITTSBURGH 
PITTSBURGH,  PA  15213 

01CY  ATTN:  N.  ZABUSKY 

DIRECTOR  OF  RESEARCH 
U.S.  NAVAL  ACADEMY 
ANNAPOLIS,  MD  21402 
02CY 
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PLEASE  DISTRIBUTE  ONE  COPY  TO  EACH  OF  THE  FOLLOWING  PEOPLE  (UNLESS  OTHERWISE 
NOTED) 


NAVAL  RESEARCH  LABORATORY 
WASHINGTON,  D.C.  20375 

DR.  H.  GURSKY  -  CODE  1)100 
DR;  P.  GOODMAN  -  CODE  4180 
DR;  P;  RODRIQUEZ  -  CODE  4706 

A.F.  GEOPHYSICS  LABORATORY 

l.g;  hanscom  field 
BEDFORD,  MA  01731 
DR.  T.  ELKINS  ' 

DR.  W.  SWIDER 
MRS.  R.  SAGALYN 
DR.  J.M.  FORBES 
DR.  T.J.  KENESHEA 
DR.  W.  BURKE 
DR.  H;  CARLSON 
DR.  J.  JASPERSE 
Dr.  F;J.  RICH 
DR.  N.  MAYNARD 

BOSTON  UNIVERSITY 
DEPARTMENT  OF  ASTRONOMY 
BOSTON,  MA  0221 5 
DR.  J.  AARONS 

CORNELL  UNIVERSITY 
ITHACA,  NY  14850 
DR.  W.E.  SWARTZ 
DR.  R.  SUDAN 
DR.  D.  FARLEY 
DR.  M.  KELLEY 

HARVARD  UNIVERSITY 
HARVARD  SQUARE 
CAMBRIDGE,  MA  02138 
DR.  M.B.  MCELROY 
DR.  R.  LINDZEN 

INSTITUTE  FOR  DEFENSE  ANALYSIS 
400  ARMY/NAVY  DRIVE 
ARLINGTON,  VA  22202 
DR.  E.  BAUER 

MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
PLASMA  FUSION  CENTER 
LIBRARY,  NW1 6-262 
CAMBRIDGE,  MA  02139 


NASA 

GODDARD  SPACE  FLIGHT  CENTER 
GREENBELT,  MD  20771 
DR.  R.F.  BENSON 
DR.  K.  MAEDA 
Dr;  S;  CURTIS 
Dr.  M.  DUBIN 

COMMANDER 

NAVAL  AIR  SYSTEMS  COMMAND 
DEPARTMENT  OF  THE  NAVY 
WASHINGTON,  D.C.  20360 
DR.  T.  CZUBA 

COMMANDER 

NAVAL  OCEAN  SYSTEMS  CENTER 
SAN  DIEGO.  CA  92152 

MR.  R.  ROSE  -  CODE  5321 

NOAA 

DIRECTOR  OF  SPACE  AND  ENVIRONMENTAL 
LABORATORY 
BOULDER,  CO  80302 
DR.  A.  GLENN  JEAN 
DR.  G.W.  ADAMS 
DR.  D.N.  ANDERSON 
DR.  K.  DAVIES 
DR.  R.  F.  DONNELLY 

OFFICE  OF  NAVAL  RESEARCH 
800  NORTH  QUINCY  STREET 
ARLINGTON,  VA  22217 
DR.  G.  JOINER 

LABORATORY  FOR  PLASMA  AND 
FUSION  ENERGIES  STUDIES 
UNIVERSITY  OF  MARYLAND 
COLLEGE  PARK,  MD  20742 
JHAN  VARY AN  HELLMAN, 

REFERENCE  LIBRARIAN 


PENNSYLVANIA  STATE  UNIVERSITY 
UNIVERSITY  PARK,  PA  16802 
DR.  J.S.  NISBET 
DR.  P.R.  ROHRBAUGH 
DR.  L;A;  CARPENTER 
DR.  M;  LEE 
DR;  R;  DIVANY 
DR.  P.  BENNETT 
DR;  F;  KLEVANS 

PRINCETON  UNIVERSITY 
PLASMA  PHYSICS  LABORATORY 
PRINCETON,  NJ  08540 
DR.  F.  PERKINS 

SCIENCE  APPLICATIONS,  INC. 

1150  PROSPECT  PLAZA 
LA  JOLLA,  CA  92037 
DR.  D.A.  HAMLIN 
DR.  L;  LINSON 
DR.  E.  FRIEMAN 

STANFORD  UNIVERSITY 
STANFORD,  CA  94305 
DR.  P.M.  BANKS 

U.S.  ARMY  ABERDEEN  RESEARCH 
AND  DEVELOPMENT  CENTER 
BALLISTIC  RESEARCH  LABORATORY 
ABERDEEN,  MD 

DR.  J.  HEIMERL 

GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AK  99701 
DR.  L.C.  LEE 

UNIVERSITY  OF  CALIFORNIA 
LOS  ALAMOS  SCIENTIFIC  LABORATORY 
J-10,  MS-664 
LOS  ALAMOS,  NM  87545 
DR.  M.  PONGRATZ 
DR.  D.  SIMONS 
DR.  G.  BARASCH 
DR.  L;  DUNCAN 
DR.  P.  BERNHARDT 
DR.  S;P.  GARY 

UNIVERSITY  OF  ILLINOIS 
DEPT.  OF  ELECTRICAL  ENGINEERING 
1406  W.  GREEN  STREET 
URBANA;  IL  61801 
DR.  ERHAN  KUDEKI 


UNIVERSITY  OF  CALIFORNIA, 

LOS  ANGELES 
405  HILLGARD  AVENUE 
LOS  ANGELES,  CA  90024 
DR.  F.V.  CORONITI 
DR;  C.  KENNEL 
DR.  A.Y.  WONG 

UNIVERSITY  OF  MARYLAND 
COLLEGE  PARK,  MD  20740 
DR.  K.  PAPADOPOULOS 
DR;  E;  OTT 

JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  MD  20810 
DR.  R.  GREENWALD 
DR;  C.  MENG 

UNIVERSITY  OF  PITTSBURGH 
PITTSBURGH,  PA  15213 
DR.  N.  ZABUSKY 
DR.  M.  BIONDI 
DR.  E.  OVERMAN 

UNIVERSITY  OF  TEXAS 
AT  DALLAS 

CENTER  FOR  SPACE  SCIENCES 
P.O.  BOX  688 
RICHARDSON,  TEXAS  75080 
DR.  R . ( HEELIS 
DR.  W;  HANSON 
DR.  J.P.  MoCLURE 

UTAH  STATE  UNIVERSITY 
4TH  AND  8TH  STREETS 
LOGAN,  UTAH  84322 
DR.  R.  HARRIS 
DR.  K.  BAKER 
DR.  R.  SC HUNK 
DR.  J.  ST. -MAURICE 

PHYSICAL  RESEARCH  LABORATORY 
PLASMA  PHYSICS  PROGRAMME 
AHMEDABAD  380  009 
INDIA 

P.J.  PATHAK,  LIBRARIAN 


